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I. INTRODUCTION 

The etymology of the word eikon traces it back to the word eikenai which is the transliteration of the word etnevai 
[1] in the Greek language meaning ‘to resemble’ [2]. In the Greek language it evolved into the word eikon which is 
the transliteration of the word clkov [1] meaning ‘image’ [2]. It was borrowed into the Latin language and later into 
the English language as eikon. In the English language it transformed into the words icon and ikon which are the 
variants of the word eikon. Notice that throughout its evolution it has meant ‘image.’ 

The main goal of this report is to understand the eikonal approximation in quantum mechanics and quantum 
field theories. Approximations play a very important role in the understanding of processes that cannot be solved 
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exactly. The Born approximation in quantum mechanics is an example of an approximation that has been extensively 
used for studying low energy processes. In quantum field theories involving coupling constants smaller than one 
we use the standard weak-coupling perturbation series which is parallel in its approach to the Born approximation 
in quantum mechanics. In the 1950’s and 1960’s when high energy physics was ascending towards its peak, it was 
realized among the high energy physicists of those times that the Born approximation is not a valid approximation for 
studying processes involving high energies. This period in fact was the golden age in the development of the eikonal 
approximation in quantum mechanics and quantum field theories. The descendents of this era took the theory of 
eikonal approximation for granted in quantum mechanics and quantum field theories. There was prolific activity in 
the application of eikonal approximation in high energy physics, especially in QCD. 

The eikonal approximation was not born in the study of quantum mechanics. It originated far back in optics. 
Light we know obeys Maxwell’s equations. In terms of Maxwell’s equations light is understood as a wave obeying a 
wave equation. But two hundred years before Maxwell wrote down his equations scientists understood reflection and 
refraction of light which was extensively studied under the branch of science called ray optics. Today in elementary 
school we learn ray optics without ever introducing Maxwell’s equations. In ray optics we assume that light travels 
in a straight line. This assumption works fine as long as the size of the obstacle is large compared to the wavelength 
of light. This is called the eikonal approximation in optics. Here we make contact with the meaning of the word 
eikon in the sense that eikon meaning ‘image’ is formed by light only in the straight line approximation. In processes 
involving diffraction we encounter the limit of the validity of the eikonal approximation in optics. We quickly switch 
over to Maxwell’s equations to study diffraction. 

Optics is described by Maxwell’s equations which can be written as a wave equation with the dispersion relation 
given by u> = kc. Similarly, quantum mechanics is described in terms of Schrodinger’s equation which is a diffusion 
equation (in imaginary time) with a dispersion relation given by u) = t^-. With this correspondence we ask the 
question, can we not have a corresponding eikonal approximation in quantum mechanics? Yes we can. The eikonal 
approximation in quantum mechanics works for processes involving the scattering of particles with large incoming 
momentum and when the scattering angle is very small. 

In the language of differential equations, the main advantage the eikonal approximation offers is that the equations 
reduce to a differential equation in a single variable. This reduction into a single variable is the result of the straight 
line approximation or the eikonal approximation which allows us to choose the straight line as a special direction. 
The early steps involved in the eikonal approximation in quantum mechanics are very closely related to the WKB 
approximation in quantum mechanics. The WKB approximation involves an expansion in terms of Planck’s constant 
h. The WKB approximation also reduces the equations into a differential equation in a single variable. But the 
complexity involved in the WKB approximation is that this variable is described by the trajectory of the particle 
which in general is complicated. The advantage of the eikonal approximation is in the classical trajectory being a 
straight line. Thus in this manner the eikonal approximation is a very stringent semi-classical limit. 

A very comprehensive collection of work on scattering theory in general with a very extensive bibliography which 
covers scattering theory in both electromagnetism and quantum mechanics is the book by Roger G. Newton [3]. A 
couple of textbooks among many which I have used are [6] and [7] . 

An extensive list of literature on the theory of eikonal approximation in quantum mechanics and quantum field 
theories developed between the years 1950 and 1970 is available in [4]. It is unanimously accepted by everyone that 
R. J. Glauber’s lecture notes on eikonal approximation [5] is the best available work on the subject. Glauber’s lecture 
notes addresses the question of the conditions for the validity of the eikonal approximations in quantum mechanics. 
In the section on eikonal approximations in this article starting out from his ideas we derive the conditions for the 
validity of the eikonal approximations more concretely. We verify our conclusions by evaluating the total scattering 
cross section for the simplest potential (an inverted finite spherical well) using the eikonal approximation and by 
comparing it with the result obtained by the partial wave expansion. 

After our study of eikonal approximation in quantum mechanics we concentrate on quantum electrodynamics. To 
study the eikonal approximation in quantum electrodynamics we need to first write the field theoretical equations in a 
suitable form. Schwinger’s formalism for quantum electrodynamics is the most suitable for this purpose. We develop 
Schwinger’s formalism in detail based on [15-18] in section V. 

In section VI we concentrate on electron-electron scattering. We introduce the quenched approximation. We outline 
the derivation of the electron-electron scattering in the eikonal approximation based on [19,20]. We point out the 
similarity of the result of electron-electron scattering obtained in the eikonal approximation to the non-relativistic 
scattering amplitude due to the Coulomb potential. 

The terminology used in this study report is the following 

f(0, <j>) = Scattering amplitude 
Oscatt. = Scattering cross section 
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cr a bs. = Absorption cross section 

C"tot. = CT scatr ^ + (T a bs. • 


II. SCATTERING IN QUANTUM MECHANICS 

A scattering process in quantum mechanics is described by the solution of the Schroedinger equation 

t ) = t ) ( 1 ) 

with the boundary conditions dictated by the requirement that the wave function ip'( f , t ) must have a component that 

j- 2 i 2 

involves an incident plane wave with energy E = , 2 K moving in the positive 2 direction and another component that 
involves a spherical outgoing wave. It should be emphasized that in scattering problems we do not require to 

go to zero at r — > oo. In fact the scattering amplitude which is the quantity of interest is contained in the r — > oo 
(asymptotic) part of ip'(r,t). 


~L v2+v ™ 


A. Formulation of the scattering problem 

For the case when the initial beam can be described by a state of definite energy we can write 

ip'(r,t) = e~k Et %jj[r) (2) 

where ip(r) satisfies the differential equation 

[V 2 + k 2 ] ip(r) = U ( r)ip(r ) (3) 

with the boundary conditions on ip(r) dictated by the boundary conditions on ?//(r, t). We have used the notation 
k 2 = and U(r) = . 

The differential equation for ip(r) in eqn. (3) can be rewritten as an integral equation [8] given by 

ip(r) = <j){r) + J d 3 r' Go(r, r')U(r')ip(r') (4) 

where 4>{r) satisfies the potential free equation 


[V 2 + fc 2 ] (j>(r) = 0 

and the Green’s function Go(r, r 1 ) is the solution to 

[V 2 + k 2 ] Go(r, f') = 5 (3 ^(r — r'). 


( 5 ) 


(6) 


The boundary conditions on <j)(f) and Go(r, f r ) are prescribed by the boundary conditions on tp(r). The integral 
equation in eqn. (4) is called the Lippmann-Schwinger equation [9] 1 . The solutions to <f>(r) and Go(r, r') are 


m = A 0 e lk ' r + B 0 e 


-ik-r 




A 


a ik\f—r\ 


B 


0 —ik\r—r\ 


( 7 ) 

(8) 


Schwinger’s view of the paper was [10]: ‘. . . I thought the importance of the paper was the variational principle . . . the so- 
called Lippmann-Schwinger scattering equation is to me conventional scattering theory written in operator notation. Nothing 
new. But that is what everybody paid attention to.’ 
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with the constraint A + B = 1. Using eqn. (7) and eqn. (8) in eqn. (4) the Lippmann-Schwinger equation takes the 
form 


= A 0 e lk ' r + B 0 e 


-ik-'i 




A 


0 ik\f-r\ 


B 


0 —ik\r—r\ 




(9) 


Imposing the boundary conditions that the wave function consists of a component that is a plane wave moving in the 
positive 2 direction and another that is an outgoing spherical wave we get 


- 4^ 


ik ' p 1 ^ d 3 r' 


ik\r-r\ 


U(r')ip{r')- 


(10) 


As we emphasized before the information related to the scattering amplitude is contained in the asymptotic region 
of the wave function. In most of the problems of interest the potential V (r) is confined to a finite volume in space 
and the detectors are very far from the region containing the potential. For these cases we can safely conclude /«r 
and thus approximate 



We can thus write 

VVwoc(r) = A 0 e ikp - -1 f d 3 r' - e lk ^~^ U (r>^) 
47 t J r 

= e^+/(M) — 

T 

where we have set Ao = 1 so that 


(ii) 


(12) 

(13) 


f(0, 4) = J d 3 r' e-^[/(r')V’(r') (14) 

can be interpreted as the scattering amplitude, where k' = k L . The illustration of the variables involved is shown in 

fig- I- 



FIG. 1. Illustration of the various variables used in the calculation. Note that r, k’ , and k are in spherical polar coordinates 
and r is in cylindrical polar coordinates. 

Often it is more useful to denote as a function of k and k' and thus write f(0,<j>) = f(k',k). Observe that 

even though the information related to f(Q , (f>) is contained in the asymptotic region of ip(r) the only contributions to 
f{9 , (f>) in eqn. (14) comes from regions where the potential is not zero. To complete the formulation of the scattering 
problem we define the scattering cross section as 

Cscatt. = j dfl | f{0,(f>) | 2 (15) 
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and the total cross section as 


&tot. — ^scatt. H” 0"abs. 


(16) 


where <r a b s . is the absorption cross section which will be explicitly defined in eqn. (27). 

To summarize this section on formulation of the scattering problem, we have concluded that a prescription to 
evaluate the scattering amplitude is to use 

m 4>) = -dj; J d 3 r' e -^[/(r'Mr') (17) 

where ip(r) is determined by solving 

[V 2 + fc 2 ] ip(r) = U(r)ip{r) (18) 

under the boundary conditions described after eqn. (1). 


B. Identities satisfied by the scattering amplitude 

Starting from the Schroedinger equation we shall derive two identities satisfied by the scattering amplitude. The 
first, the dynamical reversibility theorem was first derived by R. J. Glauber and V. Schomaker [11] in 1953. The 
second, the optical theorem, is a statement of probability conservation in quantum mechanics written as a continuity 
equation in reference to the scattering amplitude. The optical theorem was first derived by E. Feenberg [12] in 1932. 


1. Dynamical reversibility theorem 
In the Schroedinger equation if we write the wave function as 

ip'(r, t\ k ) = k) 

where 2 mE = h 2 k 2 , we can quickly derive the following identity 

ip(r; — 1’ 2 ) V 2 ip(f; ki) - ip(r; £q) V 2 0(f; -fc 2 ) = (k\ - k%)ip(r; fci)^(f; fc 2 ). 
Evaluating the above expression on a sphere of radius r — > oo for the case | k\ |=| fc 2 | we have 


dS- 


Vv — >cjo (^i ki) V lj> r — ,oo(^i^i) l^r — vqq (r, k\) V Ipr — >oo i ^2 ) 


= 0 . 


(19) 


(20) 


(21) 


Using eqn. (13) in the above expression and evaluating the surface integrals after taking the limit r — * oo we get [5] 

f(k 2 ,k 1 ) = f(-k 1 ,-k 2 ) (22) 

where we have used the notation we defined for f(9, <p) in the paragraph after eqn. (14). Thus the scattering amplitude 
for a process going from k\ to /c 2 is identical to that of a process going from — A: 2 to —k\. This is dynamical reversibility. 


2. Optical theorem 


In a very similar manner as we did in the earlier section (in this case we take the complex conjugate of the wave 
function) we arrive at the following continuity equation 

^ + V-J(r) = s(r) (23) 

where 
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(24) 


p{r) 

J(r) 

s(r) 


i m i 2 

— \ip*(r) V i/>(r) - ip{r) V ^*(f) 
2im 

| [Im ^(f)] I fp(r) I 2 • 


Integrating the continuity equation over a sphere of radius r — > oo and noting that the derivative of the p term does 
not contribute because both the states have the same energy we get 



dS- [ , 0*_ >oo (r) V 2 ^ r _ >0O (r) - tp r ^oo (r) V 2 ip*^ 00 (f f )] = lim [ d 3 r s(r). 

1 »oo J 


(25) 


Again using eqn. (13) in the above expression and evaluating the surface integrals after taking the limit r — > oo we 
get [5] 



dcp 


/»7 r 

/ sin0 dd | f(0,<j>) | 2 — lim 
Jo 


2 TO 
00 ?i 2 fc 


d 3 r [Im U(r)l | ^>(f) | 2 = — Im/(0 = 0). 

K 


Using eqn. (15) and defining the absorption scattering cross section as 


(26) 


O 777 r 

o-abs. = - lim — 2~ / d 3 r [Im V(r)] | ip(r) 
oo h k J 


(27) 


we have the optical theorem 


cr sca tt. + cr abs . = — Im/(0 = 0). (28) 


III. PARTIAL WAVE EXPANSION 

In this article we will be concerned with approximation methods for evaluating the scattering amplitude. In 
particular we will be interested in the eikonal approximation to the scattering amplitude. After we have made an 
approximation we would like to be aware of how much we have deviated from the exact result. The closest we can 
have to an exact result in a scattering problem is the result got by the method of partial wave expansion. We shall 
thus find it very useful to use the results obtained from the method of partial wave expansion as a benchmark. 

The method of partial wave expansion breaks down the initial wavefunction into an infinite sum over angular 
momentum components labeled by the quantum number l. The contribution to the scattering amplitude from each l 
is calculated separately and called the partial wave scattering amplitude. The complete scattering amplitude is then 
obtained by summing over all the partial wave scattering amplitudes. We say it is the closest we can have to an 
exact expression because the solution cannot in general be written in a closed form. Nevertheless for an incoming 
beam of a given energy characterized by the parameter k we can always define its angular momentum as kci where 
a is the scattering length defined in section III B in terms of the partial wave scattering amplitude corresponding to 
1 = 0. It turns out that for most of the potentials the contributions from / terms very much higher than ka converge 
very rapidly and thus can be neglected. Physically this is understood as, only those beams that have sufficient initial 
energy and those that pass through the scattering length of the potential are deflected. 

If the method of partial wave expansion gives us the exact solution why do we need approximation methods? 
Firstly because for almost all potentials (even for the simple Gaussian potential) we need to depend on numerical 
computational methods to solve the radial differential equation. Secondly for high energy scattering problems we need 
to sum over a large range of l. 

For potentials involving a boundary beyond which the potential is zero it sometimes becomes possible to evade the 
need for computational methods to a significant extent. For the simplest potential of this class, that of an inverted 
finite square well, the radial differential equation mentioned above can be solved exactly. Thus we shall solve the 
inverted finite square well in detail and use it for comparison with the results obtained from other approximation 
methods. 
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A. Method of partial wave expansion 


For a spherically symmetric potential of the form 

V (r) = V(r) 

we can write the solution to eqn. (18) as 

OO 

ip(r) = J^(2/ + 1 )i l Ri(r)Pi(cos9) 
1=0 

where where Ri(r) is the solution to the radial differential equation 

MKr) + 2«(r) + U _ i(l + l) _ 1 


"7~i ~ + — rr 2 - + * - - C(r) RM = o. 

ar z r dr i r z \ 

For potentials that die off faster than -y we can write the solution to eqn. (31) in the r — » oo region to be 


lim Ri(r) = Ci cos Si ( lim ji(kr)) — sin Si ( lim rji(kr) ) . 

r — too L \r — too / \r — too / . 


Substituting eqn. (32) in eqn. (30) and using the identity 


^(21 + 1 )i l ji(kr)Pi(cos 9) 


in eqn. (14) we have 


OO 

f{9 , (t>) - — = y^(2l + l)i l Pi{cos9) Uci cos Si - 1) ( lim ji{kr)\ - Q sin Si ( lim T]i(kr)) . 

V z ' L \r — too / \r— too / . 


Using 


l r 7T" 

lim ji(kr) = — cos kr — (l + 1) — 
r — *oo kr i 2. 

— 1 L— (l+l)pifcr i ,-+(l+l)p-*fcr 

" 2 kr l + J 

1 r 7T' 

lim rft(kr) = - — sin fcr — (l + 1) — 
r—> oo kr l 2 . 


_ j-(l+ 2 )gifcr _ jl e ~ikr 

2kr l 


and equating the coefficients of e ?fcr and e lkr in eqn. (34) we get 

OO 

f(0) = — 53 ( 2 1 + l)i l Ricos 9)r (l+1) [c ie +is ‘ - 1 ] 

1 oo 

0 = — V( 2Z + lJi'^Ucosfl)*- 1 -^ \Cie~ iSl - 1] . 

;=o 

The second equation above further reduces to Ci = e J<5i which when used in eqn. (37) to eliminate Ci gives 


(29) 

(30) 

(31) 

(32) 

(33) 

(34) 

(35) 

(36) 

(37) 

(38) 


= Yk ^ {21 + 1 ) P; ( C0S ^ t e2iSl - ' ( 39 ) 

lK 1=0 

This is the standard expression for the scattering amplitude given as a sum of partial waves. Si’s are defined as the 
phase shifts in Bessel functions (which are trigonometric functions in the r — > oo limit) as introduced in eqn. (32). 

Observe that the phase shifts Si ’s in eqn. (39) are in principle determined by solving the radial differential equation 
in the presence of U(r) which can be solved analytically only for special cases. In practice the phase shifts are 
evaluated by solving the radial equation numerically. For problems involving high energies a further complication 
arises because we need to sum over a sufficiently high number of partial waves. 

For the case of a inverted square well potential it is possible to solve the radial equation exactly. In the subsequent 
section we shall write down an expression for Si for the inverted square well potential. We shall find it useful to 
compare the results later when we are doing eikonal approximations. 
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B. Scattering length of a potential 


The scattering cross section is obtained by squaring eqn. (39) and integrating over all directions. Using the 
orthogonality of Legendre polynomials it turns out to be 

a OO 

/i /Jr 

CTscatt. = Tj + 1) sin 2 (40) 

o 


The contribution to the scattering cross section from the s-wave (l = 0) is 


{<7scatt.}; = 0 — 


47T 


Sill 


2 <>0- 


(41) 


The scattering length ‘a’ is defined as the radius of the sphere that will give a contribution in a classical scattering 
problem equal to the s-wave contribution above. Correspondingly we write 


Thus we get 


{c^scatt. o — 47ra 


a = — sin Jo- 
k 


(42) 

(43) 


Observe that sin (So will require the solution to the differential equation in eqn. (31) for l = 0 which normally can only 
be solved numerically. 


C. Inverted finite spherical well potential 


For a spherically symmetric potential of the form 


we can write 


V(r) 


V{r) for r < a 
0 for r > a 


Ri{r) 


i?z<(r) for r < a 

Ci [cos Siji (hr) — sin Sir]i(kr)\ for r > a 


(44) 


(45) 


where f?;<(r) is the solution to eqn. (31) for r < a. Requiring the wave function and its derivative to be continuous 
and taking the ratio of the equations got by requiring the continuity we get 


dRi<(r) 

dr 

l 

Ri<(a) 


ka [ibKM ~ tan Si-^rn(ka)_ 

[ ji(ka ) - tan Sirji(ka)] 


The above expression can be solved for tan (5/ if i?/<(r) is known at r = a. 
For an inverted square well potential defined by 


(46) 


we further have 


3 

ii 

V for r < a 

0 for r > a 

(47) 

Ri<(r) 

= Biji (ar) 

(48) 


where a = k 


1 - with 


< 1. For this case we can write the explicit expression for the phase shifts to be 


where 


sin 2 Si 


N 2 {ka , aa) 

N 2 (ka , aa) + D 2 (ka , aa) 


N(ka,aa) = aa ji-i(aa)ji(ka ) — ka ji(aa)ji-i(ka) 
D(ka,aa) = aa ji-i(aa)r/i(ka) — ka ji(aa)r]i-i(ka). 

In fig. 2 we plot Utot. versus ka ^ for the inverted spherical well potential for ka = 500. 


(49) 


(50) 



partial wave ( ka — 500) 



kd 

k E 


FIG. 2. a tot. versus ka ^ for ka = 500 using partial wave method. 


IV. EIKONAL APPROXIMATION IN QUANTUM MECHANICS 

In this section we shall convince ourselves that the eikonal approximation is valid for processes involving small angle 
scattering and very large incoming momentum. More rigorously the conditions are ^ <^1 and yj^ <C ka <C (v/e) 2 • 
In this stringent parameter zone the expression for scattering amplitude takes the form 


ne A) = - Me 


2 j / ik' -b‘ ’ 


oix(b') 


- 1 


where 


^‘-wwSZ dz ' v{9 ’ /y 


The scattering cross section takes the form 


CTscatt. = 87 ra / dtt sin" 1 

J o 

We shall derive the above formulas in this section. 


k_ 

E 


p+oo 


dz' V(b', z') 


(51) 


(52) 


(53) 


A. Formalism for the eikonal approximation 

In eqns. (17) and (18) we found the prescription for calculating the scattering amplitude to be 

/(M) = d 3 r ' e-^U^)^) 

where ip(f) is determined by solving 

[V 2 + fc 2 ] = U(r)ijj(r). 

To solve for let us write 

V>(r) = e i£ ' r >(r). 

With the above substitution and with the choice k along the z direction eqn. (55) takes the form 


(54) 


(55) 


(56) 
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2 ik— - U(b , z) 0 ( 6 , z) = — V 2 0 ( 6 , z) 
oz 

where we have used the notation r = (6, 2;). We can write the formal solution to eqn. ( 57 ) as 


0(6, 2) = 77(6, 2) — / d 2 b' / dz' G e (b,z,b',z') V' <j)(b',z ') 


where 77(6, 2) satisfies 


2i£; — U(b,z) 77(6, 2) = 0 


and G e (b, z, 6', z') satisfies 


2 ik^--U(b,z) G e (b,z,b',z') = 5 ( - 2 \b-b') S(z-z') 
oz 


with boundary conditions on 77(6, z) and G e (6 , 2, b' , 2') prescribed by the boundary conditions on 0 (r). The solutions 
to eqns. ( 59 ) and ( 60 ) are 

77(6, 2) = ■£» du U{b ’ u) ( 61 ) 

where we imposed the boundary condition 770 (6) = 77(6, 2 = —00) = 1, and 

G e (b, z, 9, z') = < 5 (2) (5 - 6') 0{z - z ') e ** /.* C/(b,u) ( 62 ) 

= 2^ ^( 5 -?) 0(2-2') 77(6,2) 77- 1 (6, 2'). ( 63 ) 

Using eqns. ( 61 ) and ( 63 ) in eqn. ( 58 ) we have 

0(&, 2) = 77(6, 2) 1 - ^ ^ ( V ^ + ^2) ^ Z ‘ ( 64 ) 

which after iteration takes the form 


0 ( 6 , z) = rj(b, 2) 1 + J dz'K{b,z',W b ,—) + J dz'K(b,z',\7 b ,—) J dz”K(b,z",\7 b ,-^j) + . . . ( 65 ) 

where the expression for K(b, 2, V&, J7) acting on an arbitrary function 3(2) is given by 


K(b,z,V b ,—)g(z)= — r] 1 (b, 2) + — j 77(6, 2) g{z). (66) 

Using the above series expansion for 0 ( 6 , 2) in eqn. ( 56 ) and plugging it in eqn. ( 54 ) we can write the scattering 
amplitude as 

m 0 ) = / (0) (0, 0 ) + / (1) { 0 , 0 ) + / (2) ( 0 , 0 ) + • . • ( 67 ) 

where 

/ (O) (fl,0) = ~^J d 2 b' £°°dz' e^-^WV) 77(6', 2') 

/ (1) ( 0 , 0 ) = / ci 2 6' [ + dz' e^-^Uftiz') r]{b',z') [ dz"K{b',z") 

4 tt 7 V-00 J-oo 

1 C /*+ 00 _ _ /»2 /*z 

/ ( 2 )( 0 , 0 ) = - — J d 2 b' j dz' e^ k ~ k) - r U(b',z') 77 (6', 2') J dz" K[b' ,z") J dz"'K{b',z'"). (68) 

We have used K(b,z,V b ,-§^) = K(b,z) for compactness. The power in the exponent *(fc — (?) • r 1 in the above 
expressions can be explicitly written as 

0 

i(k-k')-r = —ikb' sin 0 cos (0 — </>') + ikz'2 sin 2 -. ( 69 ) 

The coordinates used are pictorially described in fig. 1 . 
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B. Validity of the eikonal approximation 


Starting from the expression for in eqn. (68) and switching to dimensionless variables z —> at, b —> aw 

and V(b, z) — > Vv(b, z) where ‘a’ is the scattering length, we can write 


/(M) = ~^ {ka)2 i 


(70) 


g w {ka) + ^ g { 2\ka) + ka g^ika) + ... 

where we have suppressed the 9 and (j> dependence in jW’s. The explicit expressions for the g^°\ g ^ and g^ are 

/ /*+oo 

d 2 w' / d£ e -ikaw' sine co S ^- 4 >')+ikat' 2 s in 2 & v ^^ t ^ 

J — oo 

/ r+oo 

d 2 w f dt'\v(w',t')\ (71) 

J — OO 

g W = -- f d 2 w' f + dt' e - ikaw ’ sin0co <*-*') +ikat ' 2sin2 % v (vl\t') 2 V(w',t') 

^ J J — oo 

-jfd 2 w l f + dt' e - ikaw ' sineco <*-^+ ikat ' 2sin2 h{w',t') v(vf,t') f dt" j dt"’ \V 2 w ,v{w' ,t"')\ 

^ J J— oo */ — oo J — oo 

i n r+oo -1 /* r+oo rt 7 rt^ 

\g { A\<j d 2 w' dt' \v(w' ,t')f + - d 2 w' dt' \v(w',t')\ / dt" dt'" \V 2 w ,v(w' ,t'")\ 

= ~l [ d 2 w' [ + dt' e -ikaw' sine co S (4>-4>')+ikat' 2 sin 2 ^ ^ f dt " ,t") 2 

® J J — OO *7 — oo 

f d 2 w' f + dt' e -i‘»»'« in « c »(*-*')-H*<rt'J'i n2 l B (if )( /) j ? (uf,i') f dt" f dt'"V 2 w ,v(w' ,t'”) 

® J J — oo J — OO J — oo 

I I [ d 2 w' f dt' \v(w',t')\ f dt" \v(w',t")\ 2 

® J J — oo J — oo 


i r r+oo pt /*£ 

+ - / g? 2 w/ / dt' \v(w',t')\ / dt" / dt'"X7^,,v(w' ,t'") 

® J J — oo J — OO J — oo 


(72) 


It is easy to verify from the above expressions that for potentials that fall away faster than the contributions to 
<+)’ s are finite. For these potentials it can thus be concluded that major contribution to the scattering amplitude 
comes from the zeroth order for ^ <C 1 and ka (^) Cl. 

We shall later observe that the expression for the scattering amplitude obtained from the zeroth order eikonal 
approximation satisfies the optical theorem only in the limit fcc+ — * oo. This requirement can be achieved only if we 
additionally require ka yj g. 

Thus to summarize, the conditions under which the eikonal approximation holds are 




and 


1 , 1 

<C ka <C 


V/E 1 




(73) 


C. Zeroth order eikonal approximation 


Thus for potentials that die out sufficiently fast and for conditions | < 1 and ^g <C ka <C 
order contribution to f(9,<f>) is the most significant. From eqns. (67), (68) and (61) we can write 

r+oo 


the zeroth 


/(°) (<?,<(,) = d 2 b' J dz' e+- £, )' f '[/(6 , ,+ ?7++) 


1 

47T 


' — oo 
r+oo 


[ d 2 b' [ dz' e ~ ikb ' s i n 0cosO-0')+ifcz'2sin 2 § A g 2ife f_ cx} du U(b',u) 

J J — OO 


(74) 
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For large incoming momentum and for small scattering angles most of the contribution to the scattering cross section 
comes from angles less than or of order Thus in the zeroth order eikonal approximation we can approximate 


— ikb' sm6cos((j) — <f>’) + ikz'2 sin 2 - « — ikb'8 cos(</> — </>'). 


Using eqn. (75) in eqn. (74) we have 


(75) 


-| POO r*27T 

’"‘M-sh b ' db 'l 

We notice that the approximation involved in eqn. (75) allows us to carry out the z integral by substituting 
f-oo duU{b> , u) for a new integration variable. After evaluating the 2 integral we get 


d(fj 


/ ikb' 0 cos((p— 4> r ) 


f + °° dz'U{b', z') e ™ f-oo du U{6> ' u \ 
J — oo 


(76) 


where 


/ (0) (M) 


k 

2i ri 





d(f> 


/ g— ikb' 6 cos((p— <fr') 



X(b') 


kl_ 
2 E 



dz' V(b',z'). 


(77) 


(78) 


For the case when the potential is independent of variable such that the scattering is symmetric about the z-axis 
we can further carry out the (j)' integral and thus get 


/ (0) (<9) = - / b'db' J 0 {kb' 9) fe ix(fc,) - 1 


(79) 


where we have used the identity 


J o(t) = ^ d< t> e ~ 


(80) 


and Jo(t) = Jo( — t). For making the above expressions more illustrative we scale the integral variables with respect 
to the scattering length of the potential defined in section III B by introducing the dimensionless integral variables u 
and t defined as z = au and b = at , where a is the scattering length of the potential. We also use V to signify the 
maximum value of the function V(f). In term of these variables scattering amplitude in eqn. (79) takes the form 


where 


f(9) = a - ka t dt Jo(t, kad) 


,ika^£(t) _ i 


1 1 f + °° 

t(t) = ~2 yj du V (at, au). 


(81) 


(82) 


Using eqn. (15) we have the expression for the scattering cross section to be 


^scatt. 
7 TO 2 


2 (ka) 2 



e +ika%£(t) _ i 


e ~ika%£(t') 


1 f sin0 d9 J 0 (t ka9)Jo(t' ka9). 

- 1 Jo 


(83) 


We said earlier that ^ is a good estimate for the upper limit of the scattering angle 9. Thus in the limit | C 1 we 
can write sin 9 ~ 9 and limit our range of integration in 9 from 0 to ^ . Further introducing the integration variable 
x = ka9 we have 


^scatt. 

ira 2 


= 2 


t dt 


t' dt' 




p ka 


x dx Jo(tx)Jo(t'x). 


(84) 


Once we have the above expressions for the scattering amplitude in eqn. (81) and scattering cross section in eqn. 
(84) we would like to check if the optical theorem is satisfied in the eikonal approximation. For potentials that are 
not complex the optical theorem can be stated as 
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Using eqn. (81) in eqn. (85) we have 


Cscatt. = — Im /( 0). 

k 


(85) 


^scatt. 
7r a 1 



( 86 ) 


Comparing eqn. (84) and eqn. (86) we conclude that under the conditions |<Cl and ka (|() <1 alone the eikonal 
approximation does not satisfy the optical theorem. The physical content of the optical theorem is the statement of 
probability conservation in quantum mechanics. Can we save the eikonal approximation from violating the optical 
theorem? Yes. We recognize that in eqn. (84) if we take the limit ka — > oo and use the identity 


x dx Jo(tx)Jo(t'x) = ^ — - (87) 

we get exactly the expression required for the optical theorem to be satisfied. Thus we observe that even though 
the expression for the scattering amplitude in the eikonal approximation is derived under the conditions | <C 1 and 

ka(jj) 2 « 1 alone, we still need to put an additional condition ka ^ » 1 for it to satisfy the optical theorem. Thus 
we have the eikonal approximation valid under the conditions 



V 

E 


< 1 


and 


1 

V/E 


< ka < 


1 

(W' 


(88) 


For a given ka, the second inequality above puts an upper and lower bound on -j) given by ^ <C ^ <C -^= . This can 
be rewritten in the form 


V , — 

1 <C ka— <C Vka. 
Jb 


(89) 


D. Examples 


1 . Inverted finite spherical well potential 


An inverted finite spherical well potential is defined as 


V(r) 


V for r < a 
0 for r > a 


(90) 


Noting that b 2 + z 2 — r 2 and thus integrating 2 from 0 to \Jr 2 


b 2 we have for this potential 


£(*) = -Vi-t 2 


l f 

f(9) = a — ka t dt J 0 (t kad) 

1 Jo 



^scatt. 

7ra 2 


= 8 t dt sin 2 


V / 

ka —\/ 1 - t 2 
Jb 


(91) 


In figs. 3 and 4 we plot tr sca tt. versus ka^j and compare it with the plots got using the partial wave method. Using 
eqn. (89) for ka = 50, we have the region of validity to be 1 <C ka ^ <C 7, and for ka = 750, we have the region of 
validity to be 1 <C kaj= <C 27. These estimates for the region of validity are in agreement with the region where the 
curves obtained from the eikonal approximation fit the curve obtained from partial wave method in figs. 3 and 4. 
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C’tot. 

7ra 2 



FIG. 3. Comparison of eikonal and partial wave method ( ka = 50) for the inverted finite spherical well. 


C^tot. 

7ra 2 



FIG. 4. Comparison of eikonal and partial wave method (ka = 750) for the inverted finite spherical well. 


2. Yukawa potential 


Consider the Yukawa potential 


V(r) = -V 


c-r/a 


where V = — and r 2 = b 2 + z 2 . For this potential we get 

• u 


1 p+OO -\Zt' 2 + U 2 

1 du 


Vt 2 + U 2 Jo 


dO e 


—t cosh 6 


= K 0 (t) 


where Ko(t) is modified Bessel function of order zero. Using this in eqn. (81) we have 

1 r°° i 

f(9) = a — ka / t dt J 0 (t kaO) e ika ^ Ko{t) - 1 . 

* Jo 1 J 


(92) 


(93) 


(94) 


14 


The contribution to the above expression from Jo(t kci9) dies off very fast. Thus to a good approximation we can say 
that we get non-zero contribution to the integral for t had <C 1. Noting that in the eikonal approximation 9 < ^ and 
1 <C fca-jy <C pyjj we conclude that non zero contributions to the integral comes from t < Cl. In the limit tC 1 the 
modified Bessel function takes the form 


K 0 (t) 


In - - 7 

t ' 


(95) 


where 7 = 0.577 ... is the Euler’s constant. Further we observe that one of the terms in eqn. (94) contributes only at 
9 = 0 and thus is a delta function. Overall after substituting x = ka9 we get 


1 f°° 

f{9) = ika 2 5{9) - ; e ifc <On 2 + i„(fca 0 )- 7 ) / dx jt x \ x l ~ ika T: . 

k0 J 0 


In terms of the gamma functions we have 



dxJ n (x)x 1+2ia = 2 2ia+1 


r(i + + iot) 

r(§ -ia) 


Thus we can write 


1 m ie#o= 


1 

2 

r 1 


2 k(l) 2 

ka -g 

r 1 

( 1 + i fe 4 ) 


(96) 


(97) 


(98) 


1 2 2 

where ^ ka = j- using V = 7 - and Tik = p = mv. Observe that the above result is identical to the scattering 
amplitude due to a Coulomb potential for small scattering angles using the approximation sin0 « 9. It is worth 
observing the interesting fact that the above result is independent of ‘a’. 


V. QUANTUM ELECTRODYNAMICS 

A. Schwinger’s functional differential equations for quantum electrodynamics 

We define the vacuum to vacuum persistence amplitude for electrodynamics in the presence of external sources 
J(x) = {J^(x),ri(x),fj{x)} to be 


Z[J\ = (0+ I 0-) J . (99) 

Using Schwinger’s quantum variational principle we can write 

8Z[J\ = * (0+ | 8W[A, J] | 0-) J (100) 

where IU[>4; J\ is the action for quantum electrodynamics given in terms of the fields A(x) = 
{A^^x), F M „(x), t/)(x), tp(x)} interacting with external sources is J(x) 

W[A; J] = S[.4] + J d A x \J ,l (x)A ll (x) + fj(x)ip(x) + ip(x)ri(x)] (101) 


where 


S[A] 



- l -F^(x)(dMx) 


d„A^(x)) + ^F^(x)F /Jj „(x) + - m + A^{x)\ij}(x) 


(102) 


Gauge invariance requires us to constrain the external source J^(x) to satisfy the condition d^J^x) = 0. Using 
Schwinger’s quantum variational principle to vary Z\J\ with respect to the fields A and the external sources J we 
get 
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SZ[J\ = i6fj(x)(0 + | dW s [ ^\ J] | 0-) J + i8 V (x)(0+ \ dW } A ^ ] \ 0-) J + i 5J»(x)(0+ \ I 0_ 


i Sip(x)( 0+ | 
i SA^(x)(0+ 


5rj(x) 

8W[A,J] . 
S’ip(x) 
SW[A,J\ 


0_) 7 + i 5i/>(x){ 0 -| 


5r](x) 

(5V>(a;) 


SJ^(x) 


0_ 


..7 


O.}' 7 + * <Jf’»“'(®)(0+ | 5W } A ] J ^ \ 0-) J 


6A»{x) 1 _/ w ' ^ 1 8Fv i ’(x) 

= i 8fj(x)( 0+ | I 0-} 7 " + i <5?y(a;)(0 + | •ip( x ) | O-) 7 " + i 5J fl (x)( 0+ | A^x) | 0-} 7 " 

+ i [(i7 M 9 M ~ m)(0+ | V>(z) I 0-) J + e(0+ | 7 ^Afj,(x)ip(x) \ 0-) J + g(x)(0+ | 0_) 7 ] 

+ i (0+| ip(x) | 0_) 7r (-*7 M< a p - to) + e(0+ | %[){x)^ A^x) | O-)- 7 + fj(x)( 0+ | O-) 7 " Sijj(x) 

+ i 5A»(x) [-9 M (0+ | F^(x) | 0-) J + e(0+ | ^(x^^x) \ 0-} J + J li (x)( 0+ | O-) 7 ] 

+ * 5F^{x) 1 - [(0+ | F^{x) | 0_) 7 - {9 M (0+ | A v (x) \ 0-) J - 9,(0+ | A^x) | 0_) 7 }] . (103) 

In the above equation we have treated the variations in the fields to be c-number variations. Schwinger’s quantum 
action principle states that Z\J] is stationary with respect to variations in the dynamical parameters (fields A(x) in 
the present case). Thus we can write 

H)H|y = <»♦ i A ^> i °-> J 

H)® = (° + |^)l °-)‘ 7 

orj(x) 

(- i )§Q = (°+\t(x)\0-) J 


Srj(x ) 


0 = (i 7 M 9 M - to) ( 0 + | ip(x) | 0_) 7 + e( 0 + | 7 M A M (x)V>(a;) \ 0_) J + 7i(x)(0 + \ 0_) J 
0 = (0+ | ip(x) | 0 _)‘ 7 (-i 7 Ai< 9 M - to) + e(0+ | ip{x)Y A^x) | 0 _) 7 + 17 ( 0 ;) (0+ | O -) 7 
0 = -9^(0+ | F M ,( x) | 0 _) 7 + e(0+ | t/>(x) 'y^.ip{x) | 0 _) 7 + J^(x){ 0+ | 0 _) 7 
0 = ( 0 + I i^(*) I 0_) 7 - {9 m (0+ I A v (x) I 0_> 7 - 9,(0+ I ^(a) I 0 _) 7 }. 


Eliminating the fields in the above equations we get 


*7^9^ — m + e 7 M (— *) 


8Ji J -(x) 


^ 5 -m=’ 1{x)zlJ] 


5rj{x) 


— *7^ 9 M — to + 67^ (— t) 


5J^( x) 


= V(x) Z[J } 


g^d 2 - d^d" 


H)IS§ 


5J /J ‘(x) 


5r)(x) 


5r)(x) 


(104) 

(105) 

(106) 
(107) 


These are the functional differential equations for Z[j7"]. This was first written down in this form by Julian Schwinger 
[13,14], Since Z[j/] is the generating functional for the Green functions it has the information regarding all possible 
processes in quantum electrodynamics. Conversely a solution to the above functional differential equations solves 
quantum electrodynamics completely. 


B. A formal solution to Schwinger’s functional differential equations 


Using Schwinger’s quantum variation principle in eqn. (100) for the variation of Z\J\ with respect to e we have 

H) ?m = < 0+ 1 1 


9e 


9e 


= / d 4 x (0+ | ip(x)^A IJ ,(x)tp(x) | 0_ 




(108) 
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Using eqn. (104) to eliminate the fields we get the first order differential equation 


dZ[J\ 

de 


= i J d*x (-i) 




5rj(x) 


5J^(x) 6r)(x) 


Which can be integrated to yield the result 

Z[J\ = e ie f d4x Z e=0 [J] 


(109) 


( 110 ) 


Where Z c —q\J\ is the generating function for the case when the photon field is not coupled to the fermion field. It is 
fairly straight forward to show [19] that the solution to eqns. (105), (106) and (107) for the e = 0 case is 


/ ^ 4;C f d 4 A J l _ l (x)D t ^ (x—x')J u {x / ) J d 4 x J d 4 x' r](x)S+ (x— x')r}(x') 

where the Green functions (x — x') and S+(x — x') satisfy the differential equations 

-{g^d 2 - d^d v )D^{x - x') = g' iV 8^\ x - x') 

— (7 ^dfj, — m)S + (x — x') = <5^ ( 2 ; — x') 


( 111 ) 


( 112 ) 

(113) 


with the boundary conditions on the Green functions prescribed by the requirement |Z | 2 < 1. Using eqn. (Ill) in 
eqn. ( 110 ) we have 

Z[J\ = e ~ e / d 4 * f d 4 x , T ^j8 (4) (x * f d 4 x f d 4 x'r)(x)S + (x-x')ri(x') e % f d 4 x f d, 4 x J^ l (x)D^' {x-x')J„(x') (X14) 

For Xi, yi being Grassmann variables we have the identity 

(115) 


d 

e °«™ 


■ ( X _ e 5?' r B.(i+iB)- 1 7? e Tiln(i+i.B) 


where we have used the symbolic notation, a (vector) for a 1 and A (matrix) for Aij and Tr is the trace over the i,j 
indices. Generalizing the above result for the case of functionals dependent on Grassmann variables we have 

g J d 4 x J d 4 x'j^jM(x-x , )j^ 7 j e J d 4 x f d 4 x'fj(x)N(x—x')r](x') _ d 4 x f d 4 x'fj(x)[N(S+MN)~ 1 (x,x')]ri(x) e Trln(<5+AfiV) (H6) 

where we have used the corresponding symbolic notation to suppress the integrals, Tr is the trace over both the 
spinor and coordinate index, and <5 is the Dirac delta function. Using the above identity, the expression in eqn. (114) 
simplifies to 


Z[J] = e i -f dix / dixlr>( ~ x ' >G A x y<(- i )-s7Tr] r l( x ) e L[(-i)jj T r] e i f d 4 x J d 4 x’ (x)D^ {x-x) J„ (x 1 ) 


(117) 


where 


L[A] = — Tr In G+[x, x'\ A^] + Tr In S+(x — x') (118) 

with Tr denoting the trace over both spinor and coordinate index and the Green function G+[x, x r ; A^\ satisfies the 
differential equation 

[— (*7^9^ — to) + e 7 tJ, A fJ/ (x)\ G+[x, x'\ A = 8^(x — x') (119) 

with the boundary condition prescribed by the initial condition 

lim G+\x, x'; A^] = S+{x — x'). (120) 

e— >0 

To get prepared for the next step we write down the following identity involving the Gaussian function and an arbitrary 
function F{x) in terms of the real variable x 

F A ei xbx = ei xbx e ^ b ^F(bx). (121) 

ox 

Generalizing the above identity for functionals we have 
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2) ^ j I ^ X S d 4x ' J/ j.( x )D+ 1 '(x— x , )J u {x / ) ^ J 1 d A x J d 4 x' J /Jl (x)D^ u (x— x) J u (x') 


x { e^/^/^'sra^rC* x "> 6A*(x) F[A]| ( 122 ) 


a=Jdj 

where A = f DJ stands for A p (x) = J d 4 x'D+ I '(x — x')J u (x'). Using the above identity in eqn. (117) we have 


^ J d 4 x J d 4 x' J l _ t (x)D+ u (x— x , )J ly {x r ) 


x < e 


{■ 


if d4x f d4x> 6A*(x) D^ix-x') lA S( x ) J f d 4 x J d 4 x'ri(x)G + [x,x';A ll ]Ti(x) „L[A]' 


' a=Jdj' 


(123) 


We observe that our problem of determining Z[J] for quantum electrodynamics reduces to the evaluation of the Green 
function G+[x, x'; A p ] for an arbitrary A p (x). Any kind of approximate solution to G+[x, x'\ A p ] suggests a possible 
non-perturbative approximate solution to Z\J]. 


C. Electron-electron scattering 

We can expand Z[J] around Z[ 0] and thus write 


Z[J] = Z\ 0 ] + J d 4 x r](x) 
+ J d 4 x J p {x) 


Srj(x) 
S 


z\J\ 


J = 0 


SJ^(x) 


Z[J\ 


d 4 x 
f T- 


Srj(x) 


Z[J\ 


V{x) 


J = 0 


J = 0 


where T is the quantity of interest in the electron-electron scattering process and is given by the expression 


where we have explicitly introduced the Dirac indices for clarity and 


G AB CD ( x 1 , y 1 , x 2 , y 2 ) = \ (-*) 4 tt - 


u u z\ r\ 

Sfj A (x 2) Sfj B (x l) St] A (y2 ) Sr} A (y l) 

| e M [A] G+[xi,yi; A* 1 ] G + [x 2 , y 2 \ A p ) | 


J = 0 

A =J DJ,J=0 


(124) 


T= ~y~y j d 4 Xi J d 4 x 2 J d 4 yi J d 4 y 2 VA{x 2 )y B {x 1 )yc{y 2 )yD(yi) G AB cd(xi, Vi, x 2 , 2 / 2 ) (125) 


(126) 


where we have introduced the short form for the operator 


M 


' <5 ' 
SA 



d 4 x 


d 4 x' 


SA^^x) 


D^ix-x') 


5A l '( x ) " 


(127) 


T will be interpreted as the probability amplitude for the electron-electron scattering process. Every process in 
quantum electrodynamics can be generated out of the expansion in eqn. (124). It is due to this reason that Z\J) is 
also called the generating function. 

For the electron-electron process we shall be interested in the case when the sources of the process fj(x 1 ), fj(x 2 ), 77 ( 2 / 1 ) 
and 77 ( 2 / 2 ) are interaction free states with well defined initial or final momentums. Thus we choose 


V(xi) = e tpimi u( Pl )[^ Pll , + to] 

V(yi) = e ip ’ lVl [7 V ljU + m]u(pi) 
fj(x 2 ) = e~ V2X2 u(p 2 )[^p 2 ^ + m] 

77 ( 2 / 2 ) = e ip ' 2V2 [7 V 2(U + m]u(p' 2 ). (128) 

Accordingly T will be denoted as T(pi,p' l ,p 2 ,p' 2 ). In the next section we shall concentrate on approximate solutions 
to T{jp\ , pi , p 2 , P 2 ) • 
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VI. APPROXIMATION TO ELECTRON-ELECTRON SCATTERING 


A. Quenched or ladder approximations 


To start with we shall assume that e L ^ in eqn. (126) does not contribute. Further we assume that the operator M 
in eqn. (127) act in such a way that they connect the two fermion propagators via photon lines. In the diagrammatic 
language this amounts to omitting all those graphs that involve fermion loops and keeping only those graphs that 
have photons exchanged between the fermion lines. This is called the quenched or ladder approximation. In this 
approximation eqn. (126) reads 


Gabcd{xi, Ui,X2, 2 / 2 ) 


{e Ml2 by [G + [ Xl , yi ;A»] G + [x 2 ,y 2 ; A p }}\ 

1 J A= J DJ, J= 0 


(129) 


where we introduced M\ 2 in place of M to signify that it connects the two fermion lines. Using the identities 


e Ml2 


1+ f d\e XMl2 M V2 

Jo 


(130) 


and 


6 

SA^^z) 


G+[x 1 ,yr,A p ] 


G+[x\,z-, A p ] ejf, G+[z,yi\ A p ] 


we thus get 

T(p 1 ,p' 1 ,p 2 ,p 2 ) = -i J dX e XMl2 J d 4 Zi J d A z 2 (z\ - z 2 ) 
x J d 4 xie~ ipiXl u{pi)[^p lfl + m}G + [x 1 ,zi\A^] J 
x J d A x 2 e~ ip2X2 u{p 2 ){^p 2 ^ + m}G+{x 2 ,z 2 \A%] e^ v J 


d 4 y 1 e zp i yi G+[z 1 ,y 1 ; Al][y p p' llx + mju^) 
d 4 y 2 e ip ' 2V2 G+[z 2 , y 2 \ A p ][j p p' 2fi + m]u(p' 2 ). 


(131) 


(132) 


B. Eikonal approximation to G+[x,y; A p ] 


Next step is to evaluate G + [. x , y\ A p ] in the eikonal approximation. We defined G + [x, y; A p ] in eqn. (119) to satisfy 
the differential equation given by 


[-(i7 M d M - to) + eYA^x)] G+[x,y;A p ] = <)' (4) (:r - y). 

In the eikonal approximation we have [19-22] 

G+[x,y; A p ] = i / ds e~ zms S(x — y — s—) e ie J 0 dT ™ p A ^ X ~ T ™\ 


(133) 


(134) 


Using the above eikonal approximation to G+[x,y; A p ] in eqn. (132) it becomes possible to evaluate T(pi,p' 1 ,p 2 ,p' 2 ). 
Here we shall simply state the result without the detailed steps in the derivation. We define 


T(pi,p' 1 ,p 2 ,p' 2 ) = (2n) 4 6 i4:) (p 1 +p' 1 -p 2 - p' 2 )M(p 1 ,p' 1 ,p 2 ,p , 2 ) 
and in the eikonal approximation we have the result [19] 

e 2 (pi +p 2 ) 2 T(1 - tfj) 


M(pi,p' 1 ,p 2 ,p 2 ) |= 


2m 2 (pi — p'{) 2 T(1 + i^-) ' 


(135) 


(136) 


We notice the close similarity of the above expression to the non-relativistic scattering amplitude due to the Coulomb 
potential. This is not surprising because eikonal approximation is a stringent semi-classical limit. 
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Webster’s dictionary meaning: a person who champions the less accepted cause for the sake of argument. 
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